Deterministic macroscopic quantum superpositions of motion via quadratic 

optomechanical coupling 
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We propose a scheme to prepare macroscopic quantum superpositions of motion in optomachanical 
nano- or micromechanical oscillators quadratically coupled to an intracavity field. The nonlinear 
optomechanical coupling leads to an effective degenerate three-wave mixing interaction between the 
mechanical and cavity modes. The quantum superpositions result from the combined effects of 
the interaction and cavity dissipation. We show analytically and confirm numerically that various 
deterministic quantum superpositions can be achieved, depending on initial mechanical state. The 
effect of mechanical damping is also studied in detail via the negativity of the Wigner function. The 
present scheme can be realized in various optomechanical systems with current technology. 

PACS numbers: 



Introduction. — While various quantum superposition 
states have been realized in microscopic physical systems 
such as atoms and photons [l], 0, their realization in 
macroscopic objects remains a challenging task due to 
the increasing rate of decoherence induced by coupling to 
the environment @. Yet, this situation may soon change 
with recent rapid experimental progress in quantum op- 
tomechanics — the interface between nano or microme- 
chanical oscillators and quantum optical fields Q. In 
particular, cooling macroscopic mechanical oscillators to 
their quantum ground sta te La -lal and realizing strong op- 
tomechanical coupling f§-fll| will lay a solid foundation 
for exploring macroscopic quantum states theoretically 
and experimentally in macroscopic systems [l^ - |20| . Be- 
sides their fundamental interest, these states a re p oten- 
tial resources in various quantum technologies [2lM23l | . 

This paper analyzes theoretically a new scheme to pre- 
pare macroscopic quantum superposition states of nano- 
or micromechanical oscillators coupled quadratically to 
a light field. In contrast to existing schemes that rely 
on conditional quantum measurements (l6l - fl9| or quan- 
tum state transfer [2fj| . this scheme is deterministic. It 
relies on the combined effects of two-phonon emission 
and absorption that can be engineered with quadratic 
optomechanical coupling. Pure mechanical quantum su- 
perposition states can be generated for time scales over 
which mechanical damping remains negligible. This ap- 
proach could also be applied to other physical systems in 
which degenerate three-wave mixing can be engineered. 

Model. — We consider a "membrane-in-the-middle" 
cavity optomechanical system. The mechanical oscilla- 
tor could be e.g. a suspended dielectric membrane (as 
shown in Fig. 1) Q, an ultracold atomic cloud 24 1, or a 
trapped levitated dielectric sphere [25j | . For such geome- 
tries either linear or quadratic optomechanical coupling 
can be achieved by adjusting the equilibrium position 
of the oscillator in the standing-wave cavity field, and 
purely quadratic coupling can be obtained by positioning 
the mechanical oscillator at a node or antinode of the cav- 
ity field, as demonstrated in Refs. H, H3]- In our scheme 
quadratic coupling is exploited to achieve macroscopic 



FIG. 1: Possible membrane-in-the-middle setup for the 
steady-state generation of macroscopic superpositions, with 
a pure quadratically coupled to two classical driving fields to 
generate an effective two-phonon three-wave mixing interac- 
tion. 



quantum superpositions of motion, while linear coupling 
with a weak probe field could be used to characterize the 
motional superpositions. 

We assume that a cavity mode of resonance frequency 
uj c is driven by two laser fields of frequencies uij and am- 
plitudes £j (J = 1,2). In the rotating frame with respect 
to the frequency u)\ , the Hamiltonian of the optomechan- 
ical system reads 



H a /h- 



5 C A*A 
(£i + £ 2 e 



(1) 



where A and B (A^ and W) are the annihilation (cre- 
ation) operators of the cavity mode and the mechanical 
mode of frequency ui m , respectively, the detunings are 
S c = uj c — u)i and 5i2 = 0J2 — Wi, and go is the single- 
photon quadratic coupling frequency. The last two terms 
in Eq. ([T]) describe the pumping of the cavity by the ex- 
ternal lasers. 

For sufficiently strong classical driving fields we can 
decompose the field operators A and B as the sum of 
their expectation values and small quantum fluctuations, 
A = a + a and B = (3 + b. The classical amplitudes 
a = (A) and j3 = (B) are determined by the equations 

a = -{k c + iS c )a - ig Q a{f3 + /3*) 2 - + &e-* iat ), 

(2a) 

(2b) 



= -(7m + ium)P - 2*.goM 2 (/3 + /?*). 
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where k is the cavity dissipation rate and "f m the mechan- 
ical damping rate. For 3> |£T 2 1 this latter amplitude 
can be neglected in Eq. (|2af and the steady-state cavity 
amplitude is approximately a s ~ £\/{in c — A c ), where 
A c = S c - ffo (A, + &*) 2 - That amplitude can be taken 
to be real by a proper choice of the phase of £\. The 
steady-state amplitude j3 s is determined by 



(4 3oWm |a s | 2 - < - 7 2J( A + £*) = 0. 



(3) 



Unless the first term in parenthesis is equal to zero, and 
in particular for go < 0, a situation corresponding to the 
oscillator at a minimum of the standing-wave intracavity 
field, this gives f3 s = 0. 

In the following we consider the specific situation 
where the cavity is resonantly driven by the weak laser 
field and by the second red-detuned sideband of the 
strong laser field, 



S c = uj c — lj\ = 2w„ 
uj 2 = CJ C . 



(4a) 
(4b) 



Substituting the expressions for the operators A and B 
into the Hamiltonian (TTJ, considering a mechanical fre- 
quency U) m 3> g = goa s , and neglecting the rapidly oscil- 
lating terms, the resonant optomechanical Hamiltonian 
becomes then approximately 



H = KgaP 2 + H£ 2 a + h.c. 



(5) 



where g can be controlled by adjusting the laser ampli- 
tude \£i\. This Hamiltonian can be interpreted as de- 
scribing the degenerate three-wave mixing of optical and 
acoustic waves, whereby a photon is absorbed from (or 
emitted into) the pump laser of classical amplitude £i 
and two phonons are simultaneously created (or anni- 
hilated). Due to the correlated nature of two-phonon 
emission, quantum features such as phononic squeezed 
states can be achieved in such a process [2q [. Here we 
take advantage instead of this interaction to explore the 
generation of macroscopic quantum superpositions of me- 
chanical motion. 

When taking into account cavity dissipation and me- 
chanical damping, the density operator p of the full op- 
tomechanical system is governed by the master equation 



p(t) = ±-[H,p] 

in 



£ K p + C lm p, 



(6) 



with 



C K p = n{2apa) — a) hp — pa) a), 
£-y m P = 7m ("th + l)(2bpb f - b ] bp - p 
+ j m n th (2tfpb - Wp - pbV) 



where CkP accounts for the optical cavity losses, C lm p 
describes the mechanical damping in its thermal en- 
vironment, and the mean thermal phonon number is 
n th = ( e ^n,/feBT _ w ith k B the Boltzmann con- 

stant and T the temperature. Equation ([5]) is our starting 



point to investigate the generation of quantum superpo- 
sition states of the mechanical oscillator. 

Dissipation-induced superpositions. — We proceed by 
first neglecting mechanical damping, 7,„ = 0, in which 
case it is not difficult to find that with the aid of cavity 
dissipation the full system evolves toward the dark state 



IVm) = |0 o )|Vfe,. 



(7) 



where the mechanical steady state \ipb,s) satisfies the 
equation 



{b 2 +£ 2 /g)\^ s ) =0, 



(8) 



or equivalently b \ip. 



'& )S ) = £ \ipb,s), with £ = -£ 2 /g- It 
is readily verified that the state ([7]) then satisfies the dark 
state conditions H\tpd) = and £ K (|?/'d)(V'd|) = 0. 

By expressing \ipb s) in the Fock basis as \ipb,s) — 



Er 



\rib), it follows from Eq. ([8]) that 



Cn b +2 = £oC n J\/(n b + l)(n b + 2). 



(9) 



That the probability amplitudes c nb with even (or odd) 
rib are only related to those with even (or odd) phonon 
numbers in the absence of mechanical dissipation follows 
from the two-phonon nature of the interaction. For a 
mechanical mode initially in an even Fock state we then 
find from Eq. © that 



1 



b,s 



V / cosh]£ ] 
AC 1/2 (l/3) 



E 



C 



-/?», 



\2n b ) 



(10) 



which is an even coherent state, with j3 = \/£~o and the 
normalization Af e = 2[1 + exp(— 2|/3| 2 )]. Likewise, when 
starting from an initial odd Fock state, the mechanical 
mode decays to the odd coherent state 
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b,s 1 
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c 



Vcosh\£ \ n % y/(2n b + l)\ 

Af- 1/2 m- 1 -/*», 



\2n b + 1) 



(11) 



with No = 2[1 - exp(-2|/?| 2 )]. As schematically illus- 
trated in Fig. [2l these steady-state superpositions result 
from the combined effects of cavity driving and dissi- 
pation and two-phonon emission and absorption. Such 
superpositions was proposed in a trapped cold ion via 
qubit-phonon coupling [27( . 

It is generally much easier to prepare a mechanical vac- 
uum or coherent state than a Fock state with an odd or 
even initial phonon number. In that case the mechanical 
oscillator will evolve toward the mixed state 

Pb.s = Poo\tpb,s)oo(4>b,s\ + Pu\ipb,s)n{if>b,s\ 

+ (Poil^sMVvl+k-c.), (12) 

where P o = En b ( 2n b\pb(0)\2n b ) , P n = En b ( 2n b + 
l\pb(0)\2n b + 1), and p b {0) is the initial state of the me- 
chanical oscillator. 
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FIG. 2: Schematic of the processes involved in the realiza- 
tion of the dark state \ipb,s)o in case the initial mechanical 
state contain only even phonon numbers. The straight verti- 
cal arrows describe the optical driving of the resonator, the 
wavy vertical arrows account for optical dissipation, and the 
horizontal arrows show the two-phonon transitions in the me- 
chanical oscillator. 



The purity of that state is characterized by the cross- 
term probability Poi, which depends in general on the 
initial coherence of the mechanics as well as on the sys- 
tem parameters. Assuming that £2=0 we can obtain an 
analytical estimate of its value in the bad cavity limit, 
in which case one can eliminate adiabatically the cavity 
mode, giving the reduced master equation for the oscil- 
lator 



dt 



Pb(t) = 72(2& 2 M t2 - bH 2 Pb - Pb bH 2 ). (13) 



This equation describes an effective dissipation mecha- 
nism via two-phonon absorption at the rate 72 = g 2 / K. 
We have from Eq. (fl"2"j) that the initial state p b (0) = 
I A)) (A> I, with /3q = |A)| exp(i# ), will decay to the steady 
state 



1.0 
Pb 

0.5 



0.0 



■1 


g=0.lK 

g=0.2 K 
g=0.5K 


: 


(a) 




-g=0.1 K 
g=o.2 K 

g=0.5K 



(b) 



20 40 60 80 100 20 40 60 80 100 



Kt 



Kt 




FIG. 3: Time evolution of (a) the purity Pi, and (b) the fi- 
delity Fo = \/o{ipb.s\pb\ipb,s)o for ft — VE and initial vacua of 
both the cavity and the mechanical modes. Here the pumping 
amplitude £2 = 5g and the mechanical damping 7 m = 0. (c) 
Wigner function Wb of the mechanical state in the long-time 
regime. Time in units of 1/k. 
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Pb,s = Poo|0 6 )(0 b | + Pii|l 6 )(lfe| + (Poi|0 h )(l b | +h.c), 

(14) 

where P 00 = 1 - P u = [1 + cxp(-2|/3 | 2 )]/2 and 

Pai = |A|e~ l/?o|2_4e °/o(|A| 2 ), where In(x) is the mod- 
ified Bessel function of the first kind [2{|. For the rel- 
atively large initial coherent amplitude |/?o| = 1-7, the 
probability amplitudes Poo — Pn — 0.5 and |Poi| ~ 0.42. 
That is, the state reached by the oscillator is approxi- 
mately the quantum superposition 

\^ b ,s)2^(\Q b )+e t9o \l b ))/V2. (15) 

Numerical results. — To confirm these approximate an- 
alytical results and evaluate the effects of thermal fluc- 
tuations, we solved numerically the full master equation 
(|6]), expanding the density matrix p on the Fock basis 
{\n ai n b )} as 

P = ^2 Pma,m b ,n a ,n b \n a ,nb)(m a ,m b \. (16) 



FIG. 4: (a) Dependence of the purity Pb, the probability am- 
plitudes Poo and Pn, and the coherence |Poi| of the steady 
mechanical states on the ratio g/n for initial cavity vacuum 
and mechanical coherent state |/?o) with \(3q\ = 1.7. Here 
the pumping amplitude £2 — and the mechanical damping 
7m = 0. (b) Wigner function Wb of the steady mechanical 
state for g = 0.1k. 

and truncating the system to a finite Fock space large 
enough to avoid boundary issues. We used the Wigner 
function HU 

W b(0 = i[[ d 2 V TY(e^- r '' h p b )e^'-^\ (17) 

where p b = Tr a (p) is the reduced mechanical density ma- 
trix, to reveal the quantum features of the state of the 
oscillator: negative Wigner functions are indicative of 
highly nonclassical superposition states. In addition, the 
purity of the mechanical states is characterized by the 
quantity P b = Tr(p 2 ). 

Fig.s 3 (a) and 3 (b) show the time evolution of the 
purity P b and of the fidelity P = y/o(^b,s\Pb\^b, s )o of 
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FIG. 5: (a) and (b) Time evolution of the minimal negativity 
of the Wigner function for cavity and mechanical modes 
initially in a vacuum. Here the pumping amplitude £2 = 5g 
and coupling g = 0.1k. (c) and (d) are the same as (a) and 
(b), respectively, except that the mechanical mode is initially 
in the coherent state |/3o) with /3o = 1.7 and £2 = 0. 

the mechanical state compared to the "ideal" dark state 
| , 0b,s)o for several values of the coupling strength g and 
in the absence of mechanical damping, 7 m = 0. Both 
of the cavity and mechanical modes start from vacua 
and the pumping amplitude £2 = 5g. As expected, the 
mechanical oscillator evolves asymptotically from an ini- 
tial vacuum to the pure even coherent state I'i/'b.s)!)! the 
steady-state superposition exhibiting highly nonclassical 
features demonstrated by the negativity of its Wigner 
function Wb, see Fig. 3 (c). 

The steady-state populations Poo and Pn, the coher- 
ence |Poi I j and the purity p, of the mechanical state when 
the mechanical oscillator starts from a coherent state \(3q) 
with \/3o\ = 1.7 and £2 = are shown in Fig. 4 (a). 
We find Poo — P11 — 0.5, independently of the system 
parameters. However, |Poi| decreases, from the maxi- 
mum value about 0.42 for g = 0.1k, as the coupling g 
is increased and the system departs from the bad-cavity 
limit. Likewise, the purity increases with decreasing g/ n, 
in good agreement with the analytical results. 

The effect of the mechanical damping j m on the me- 
chanical superposition is illustrated in Fig. 5. The dy- 
namics is now largely determined by two characteristic 
times: the first one is the inverse optomechanical fre- 
quency g , and the second the much longer time 7" 1 
required for the mechanics to reach a steady state. As 
would be expected the negativity of the Wigner function 
disappears over that longer time scale. Still, for mechani- 
cal oscillators with high quality factor Q = u! m /-f m , such 
as e.g. the system reported in Ref. 0, with uj m /2n = 
10 6 Hz, 7 m /27r = 0.1 Hz, and k/2tt = 10 5 Hz, highly 
nonclassical quantum superpositions of high purity can 
be still achieved in the short-time regime for tempera- 



tures as high as T ~ 10 mK (with a corresponding ther- 
mal phonon number of n t h = 200), as illustrated in Fig. 5 
(b) and (d). 

Discussion and conclusions. — In summary we have 
proposed a deterministic scheme to generate macroscopic 
quantum superposition states of a micromcchanical oscil- 
lators by quadratic cavity optomcchanics. These quan- 
tum superpositions result from the combined effects of 
optical damping and the degenerate three-wave mixing 
interaction between optical and acoustic waves resulting 
from quadratic optomechanical coupling. We showed an- 
alytically and verified numerically that the mechanical 
system can approach asymptotically a variety of quantum 
superpositions, depending on the initial mechanical state 
and pumping condition. Pure even and odd coherent 
superposition states can be obtained provided that the 
mechanical damping is neglected. It is also revealed that 
in the presence of the mechanical damping, the highly 
nonclassical superpositions can still be achieved in the 
transient regime for the realistic system parameters. 

This proposed scheme can be realized in a number of 
state-of-the-art optomechanical systems. It could also be 
used to produce quantum superpositions in other sys- 
tems where degenerate three-wave mixing can be engi- 
neered, such as microwave photons in circuit-QED con- 
figurations [30j . 

The experimental verification and characterization of 
these quantum superposition states could be carried out 
by exciting another cavity mode linearly coupled to the 
mechanical oscillator with a red-tuned probe laser. In 
the weak coupling regime, quantum-state transfer be- 
tween the cavity and the mechanical modes can then take 
place via a beam-splitter interaction (3lj . Once the me- 
chanical state is completely mapped onto the cavity field, 
one can reconstruct its Wigner function by quantum ho- 
modyne tomography via detecting the output field, see 
e.g. Ref. [32| . Future work will consider non-degenerate 
three- wave mixing in a cavity with two mechanical oscil- 
lators, with the goal to pre par e non-Gaussian entangled 
states such as NOON states |33l |. 
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Note added. — After completion of this work, we be- 
came aware of the preprint papers, arXiv:1212.4795 by 
M. J. Evcritt, T. P. Spille r, G. J. Milburn R. D. Wil- 
son, and A. M. Zagoskin; larXiv:1302.1707l by A. Vojc, 
A. Croy, and A. Isacsson, on the generation of mechani- 
cal states by two-phonon process. 
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